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INTRODUCTION 


Treatment  of  radiative  transfer  for  an  isolated  Lorentz  line  in  the 

(12) 

Lindquist -Simmons  approximation'  ' ' requires  the  evaluation  of  the 
derivative  function  y(x,  p)  defined  by 


y(x,  p)  = 


- f 

* Jo 


( 2x  ( dz 

'*'’177771777* 


or,  with  the  change  of  variable  tan(0/2)  = pz,  by 


y(x,  p) 


It  f 
* Jo 


-X(1+COS0 ) 


(p2  + 1)  + (p2  - 1 ) CO80 


d0  . 


(la) 


(lb) 


This  same  function  is  used  in  the  resonance  absorption  band  model  for- 
(3) 

mulation  of  Cobb.  ' In  both  applications,  x is  a measure  of  optical  depth, 


and  p is  a ratio  of  line  widths.  The  range  of  both  x and  p is  zero  to  infinity. 

Prior  work  on  the  evaluation  of  y(x,  p)  includes  the  rational  approx- 
imation of  Lindquist  and  Simmons^  ^ (accuracy  not  stated),  the  tabulation 
(2)  5 

of  Young'  ' (accuracy  of  one  part  in  10  ),  and  the  series  expansions  of 
Cobb^  (accuracy  of  < 1%  for  0.  01  < x < 10^,  10  2 S p < 10^,  and 

_3 

10  < y < 1).  The  present  work  considers  the  efficient  calculation  of 

y(x,  p)  to  a relative  accuracy  of  SO.  01%  for  the  entire  positive  quandrant  of 
the  xp  plane. 


The  general  features  of  the  function  are  shown  in  Fig.  1.  The  function 
limits  are 


y <x,?) 


y(x.  p)- 


l 

p/V2nx 

-2x 


e‘Xlo(x) 


x 

P 

P 

P- 


0 

1 


(2) 


where  Iq(x)  i*  the  modified  Bessel  function  of  order  zero. 

The  dramatic  difference  in  basic  mathematical  form  that  the  function 
assumes  in  different  regions  of  the  xp  plane  makes  it  impractical  to  use 
only  one  method  of  calculation.  Six  approximation  expansions  are  used  in 
the  present  method.  These  approximations  are  considered  in  the  following 
section,  and  the  regions  in  which  they  apply  are  indicated  in  Fig.  2. 
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METHODS  OF  COMPUTATION 


A.  EXPANSION  FOR  SMALL,  x (REGION  1) 

For  x < 1,  an  expansion  of  the  exponential  factor  of  equation  (la)  is 
made,  and  y(x,  p)  is  written  as 


y(x,p)=£M)ni!r-En<P> 

n = 0 


(3) 


whe  re 


E (p)  = - f ^ — =-j— 

n * •'O  (1  + z 2)  (1  + p2z2)n 


(4) 


Explicit  evaluation  of  E^fp)  by  the  method  of  residues  for  the  first  few 


values  of  n leads  by  induction  to  the  convenient  recurrence  relation 


C 2n  - 1 


En<«*>=< 


2n  n-1 


E .<P> 


P = 1 


(5a) 


* E„-1<P>J 

P * 1 


(5b) 


with 


EQ(p)  = 1 (for  all  p), 


r 


li 


and 


'W  WU'kJ'  "1FW 


1 


The  explicit  solutions  for  the  first  five  E^(p)  functions  are 


! ! 


1 

it 


E0(p)  - 1. 

Ej(p)  = l/(p  + 1). 

E,(p)  = (p  + 2)/2(p  + 1)Z,  (7) 

E3(p)  = (3pZ  + dp  + 8)/8(  1 + p)3, 

E_j(p)  = (5p3  + 20p"“  + 29p  + 16)/16(1  + p)*^. 


Additional  properties  of  E^(p)  that  can  be  derived  are 


En<P> 


v 


"n  f 1 
L 1 / P V "n 


P 

P 

P 

P 


(8) 


In  practice,  y(x,  p)  is  computed  from  equation  (3)  and  the  recurrence  relations 
equations  (5)  and  (6).  Since  equation  (3)  is  an  alternating  series  of  decreas- 
ing terms,  the  accuracy  of  the  approximation  can  be  checked  as  each  term 
is  added  to  the  series  by  testing  the  absolute  value  of  the  term  added.  The 
slowest  convergence  of  equation  (3)  occurs  for  p = 0,  in  which  case  terms 
through  n - 11  must  be  retained  for  a relative  accuracy  of  0.01%.  The 
forward  recurrence  relation  equation  (5b)  is  unstable  for  p near  unity.  For 
11  - p|  S 0.044,  this  stability  sets  in  below  n - 1 1 . However,  for  this 
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condition of  p and  for  x < 1,  y(x,  p)  may  be  computed  to  the  desired  accuracy 
with  only  the  first  two  terms  (i.e. , n = 0 and  1)  of  the  Bessel  function  expan- 
sion considered  in  the  following  section. 

B.  EXPANSION  FOR  SMALL  p (REGION  2) 

For  this  and  the  following  three  approximations,  the  starting  point 

(3) 

is  Cobb's  expansion  of  y(x,  p)  as  a series  of  modified  Bessel  functions. 

The  expansion  is 


00 

y(x,  p)  = e'X  IQ(x)  +2^  «nI  (x) 

n= 1 


where  a = (p  - 1 )/(p  + 1 ).  For  small  p,  a — - 1,  and  we  expand  y(x,  p)  in  a 

Taylor  series  about  a = -1.  At  a = -1,  the  series  converges^  to  give 
. 2x 

y(x,  0)  - e . The  first  derivative  of  equation  (9)  evaluated  at  a - -1  is 


2 

of  = -i  n=  1 


(4) 

Use  of  the  recurrence  relation1  ' 


V*'-Tn h».lW  - W-W 

in  equation  (10)  and  shifts  of  summation  indices  such  that  the  sums  over 

I . and  I , , are  transformed  into  sums  over  I yield 
n - 1 n+i  n * 

00  00 

ail  = *''*2  -£  = «'X(x>  - iiwi- 

la  = -1  n=0  n=2 
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Higher -order  derivatives  can  be  obtained  by  the  same  procedure.  The 
resulting  expansion  to  third  order  in  p for  y(x,  p)  is 


y (X , P ) 


(i: 


whe  re 


A{  xe'x(i0(x)  - Ij(x)J, 


and 


A3  = -|e'X[<2x  - 3)Iq(x)  - (2x  - 2)1^)]  - |e‘2x. 


I he  slowest  rate  of  convergence  of  equation  (13)  occurs  for  large  x. 
In  this  case,  the  coefficients  approach 


A ~ 
n 


1 


\/2itx 


(H 


and  successive  coefficients  decrease  by  the  factor  two.  For  smaller  x, 
the  decrease  is  faster.  Thus,  a worst  case  estimate  of  the  accuracy  of 
equation  (13)  can  be  obtained  by  substituting  equation  (14)  into  (13)  to  obtain 


■ 2x 


+ a + a + a 


1. 


y ■ e 


where  a = p / < 1 + p).  If  we  retain  only  terms  through  a in  the  sum,  the 
absolute  error  of  truncation  is 


e 


The  largest  a that  can  be  tolerated  for  a maximum  relative  error  of  10 
(in  the  sum)  is  given  by 


(a  + o?  + a? ) £ 1 0^  7— 

for  which  a S 0.  047  (p  < 0.  045)  is  the  solution.  Similar  analyses  show  that 

3 2 

the  o term  may  be  neglected  for  p < 0.  01,  and  the  a term  neglected  when 

p < 0.  0001. 

Computation  of  the  Iq  and  1^  Bessel  functions  required  to  evaluate 
Aj,  A^,  and  A3  of  equation  (13)  are  handled  with  rational  approximations 
that  provide  sufficient  accuracy  for  1 < x 20.  Beyond  x a:  20,  these  rational 
approximations  are  not  accurate  enough  to  yield  the  difference  Iq  - Ij 
sufficient  accuracy.  For  x £20,  use  is  made  of  the  asymptotic  expansion 

y = e~2X  +1  + 0.2bbazy  (15) 

2 

Up  to  the  q term,  this  expansion  is  obtained  from  equation  (13)  and  the 
asymptotic  expansions'  for  IQ  and  I The  coefficient  0.266  on  the  a 
term  was  determined  empirically  so  as  to  yield  0.01%  accuracy  for  all 
x 2 20  and  p < 0. 05. 
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B3  = |e‘X[(2x  + 3)I0(x)  + (2x  + 2)Ij(x)]  - |. 

| i 

As  for  the  small  p expansion,  the  worst  case  of  convergence  of 
equation  (16)  occurs  for  large  x.  In  this  case,  the  first  four  coefficients 
approach 


T" 


and 


U4  ~ 32  *2x*  ' 


Substitution  of  these  coefficients  into  equation  (16)  yields 


1 - 


2a 


1/2 


V n 


_i_  03/2 

3VW  a + 2 


where  a = 2x/(p  + 1)  . Since  the  series  alternates  in  sign,  we  can  be 
assured  that  the  maximum  absolute  error  incurred  by  truncating  equation 
(lb)  at  the  cubed  term  is  less  than  or  equal  to  a2/2.  The  condition  for 
relative  error  less  than  0.  01%  is 


a2  < 4[ 

T 1°  ' 


2a 1/2  + 8a3/2' 

VW  + ° ' \Tn 


from  which  we  obtain  a < 0.  01328 


or 


p > 12. 3X1/2  - 1. 


In  addition,  the  third;  second;  and  first-order  terms  of  equation  (16)  may 
be  neglected,  respectively,  when 


p > 28.2x1/2  - 1, 
p 2 142x1/2  - 1, 


and 


p 2 16000x1/2  - 1, 
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[ 


As  in  the  small  p approximation,  the  functions  Ig  and  Ij  required 

to  compute  the  B coefficients  are  obtained  from  accurate  rational 
(4) 

approximations. 

D.  EXPANSION  FOR  LARGE  x (REGION  4) 

Substitution  of  the  asymptotic  expansions'4^  for  I^fx)  into  equation  (9) 
and  a collection  of  terms  on  like  powers  of  x yields 


y(x,  p) 


\/2  TTX 


G(P). 

2'.(8x)<: 


(17) 


whr  re 


3 

F (p ) = fl  - 9*  - 9aZ  + or3], 

G(p)  = ) 3(3  - 25a  + 150a2  + 150c*3  - 25a4  + 3a5], 


and  a - (p  - l)/(p  + 1).  For  p ^ 1.225,  F(p)  is  negative  and  decreases 

3 

monatonically  to  F - -2p  as  p -*  <».  Similarly,  G(p)  is  positive  and 
increases  monatonically  to  G = 24p  as  p -•  °°.  The  next  higher  coefficient 
approaches  H = -720p<  for  large  p.  Thus,  the  worst  case  of  convergence 
for  p > 1. 225  appears  to  the  alternating  series 


s/Zttx 


r . / 2\  3 / 2\2  .,  / 2V3 

I 2 y2x  J 4 ^2x  j 8 y2x  J 


The  absolute  error  of  truncation  at  the  x-2  term  is  less  than  1 5(p2 /2x)3/8, 
and  the  condition  for  a relative  accuracy  of  less  than  0.01%  is 
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.c  / 2\  <1  I / 2\ 

T(fe)  5 '»  1 - j(fe) 


which  yields  the  condition  p /2x  £ 0.  0375  or 


13.  3p 


Comparisons  of  the  results  of  equation  (17)  with  exact  calculations' 

2 22 

confirm  this  bound  by  predicting  the  less  stringent  bound  x 2 7.25p  ’ "for 
20  £ x £ 104.  In  addition,  these  comparisons  show  that  x must  be  greater 
than  ~20  in  order  to  obtain  the  desired  accuracy  of  0.  01%.  This  result  is 
consistent  with  equation  (18)  and  the  lower  bound  p = 1.225  set  on  p for  this 
analysis . 

For  p<  1.22  5,  the  preceding  error  analysis  does  not  apply  because 
the  series  is  not  an  alternating  series.  Here,  comparisons  with  exact 
calculations  confirm  that  equation  (17)  is  accurate  to  the  desired  accuracy 
for  x £ 20  and  p down  to  at  least  the  upper  bound  (p  = 0.  045)  of  Region  2. 

E.  EXPANSION  FOR  LARGE  x AND  p (REGION  5) 

The  large  x expansion  equation  (17)  is  restricted  by  the  condition 
1/2 

p < 0.274x  ' , whereas  the  large  p expansion  equation  (16)  is  restricted 
by  p a 12.  3x*/2_  j.  The  following  method  provides  a means  of  computation 
in  this  excluded  region.  In  equation  (la),  we  make  the  transformation 


ITT  2 2 

1 + p z 


to  obtain 


y(x,  p)  = 


L f 

Jo  r — 


T - u 
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whe  re 


f(u) 


1 ♦ u < pfc  - 1) 


When  x is  large  and  p > 1,  f(u)  peaks  sharply  at  u 0,  and  most  of  the 
contribution  to  equation  (19)  will  come  from  this  region.  This  phenomenon 
suggests  that  not  much  error  will  result  if  we  set  the  upper  limit  of  integra- 
tion to  infinity  and  expand  the  square  root  term  of  equation  (19)  in  a power 
series.  This  latter  expansion  is 


I V'*  o 

} B u 

— — ^ n 

si  l - U2  n 0 


The  expansion  coefficient  B is 

n 


n 22n(n'.)2 


which,  for  computational  purposes,  can  be  written  in  the  recurrence  form 


B = (2n / ■ )b  . 
n \ 2n  / n - 1 


B0-  1. 


Substitution  of  equation  (20)  into  (19),  an  interchange  of  the  orders  of 
integration  and  summation,  and  setting  the  upper  limit  of  integration  to 
infinity  yields 


With  the  transformation  z 


A recurrence  relation  can  be  obtained  for  F (or)  as  follows.  Differentiate 

n 

equation  (23)  with  respect  to  a to  obtain 


From  the  solution 


Fn(o)  = e [1  - erf (0)1 


equation  (24)  can  be  used  to  evaluate  the  first  few  F^(a)  functions,  from 
which,  by  induction,  we  can  obtain  either  the  explicit  solution 


Fn(o)  = (-1)‘ 


Fo(a| 


2 

7* 


E(-n 


i (2i 


i-1 


l)(2i  - 3)  • 

72TTT7 


2 or 


(2i  - 1) 


or  the  recurrence  relations 


and 


rn<°>*Cn(*>-Fn-l<a> 


(26) 


Cn^  = 7-T-3cn-l^ 

2o 


with 


C . (o)  . 

I VITO 

In  application,  y(x,  p)  is  computed  from  equation  (22)  with  13^  evaluated  by 
equation  (21)  and  F (a ) by  equations  (2 S ) and  (2b). 

The  accuracy  of  this  expansion  was  determined  by  comparison  with 
exact  calculations.  In  the  region  x 2 5 and  p » \rZ,  the  approximation 
yields  the  desired  accuracy  with  only  a few  terms  in  the  expansion.  Near 
x *>,  p v2,  about  ten  terms  are  required,  whereas  for  x 2 10  or 
p 2 2 \ 10^  (and  x 2 1 ),  only  the  first  term  of  equation  (22)  is  required,  i.  e. , 

2 

y ~ ea  f 1 - erfar]. 


-18- 
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In  practice,  the  sum  is  terminated  when  the  relative  term  contribution  to 
the  accumulated  sum  is  less  than  10  ^ . 

Computation  of  the  function  Fn(a)  defined  by  equation  (25)  is  made  to 

-6  ^ 

a relative  accuracy  of  5 X 1 0 by  using  the  approximation  of  Matta  and 
Reichel.^ 

F.  INTERMEDIATE  x AND  p (REGION  6) 

For  intermediate  values  of  x and  p,  none  of  the  methods  considered 
to  now  can  yield  the  accuracy  desired  in  y(x,  p).  For  this  excluded  region 
(Region  6),  we  resort  to  the  Bessel  function  expansion  equation  (9) 

y(x,  p)  = e 


!0(x) 


(27) 


Although  this  expression  (with  N -*  »)  is  an  exact  solution  of  equation  (1), 
its  use  for  efficient  computation  is  limited,  especially  for  large  x,  unless 
p is  very  close  to  unity.  For  large  x,  the  functions  In(x)  fall  off  very  slowly 
with  n,  and,  consequently,  a very  large  number  of  terms  must  be  included 
in  the  summation.  The  number  of  terms  (N)  that  must  be  included  within 
Region  6 for  a relative  accuracy  of  0.01%  are  tabulated  in  Table  i.  Even 
for  x only  as  large  as  20,  as  many  as  27  terms  must  be  included  for  the 
desired  accuracy. 

In  practice,  the  number  of  terms  used  in  the  expansion  is  determined 
by  the  worst  case  of  p ~ 0.  05,  for  which 


N a 8 x 


0.  393 


The  Bessel  functions  I (x)  are  computed  with  a backward  recurrence 
17 ) n 

algorithm'  ' initiated  at  order 


M = 14x 


0.  307 


to  ensure  10"*  relative  accuracy  in  In(x)  through  n = N. 
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LABORATORY  OPERATIONS 


Th»  Laboratory  Operations  of  The  Aerospace  Corporation  is  conducting 
experimental  and  theoretical  investigations  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  military  concepts  and  systems.  Ver- 
satility and  flexibility  have  been  developed  to  a high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  spat  e and  missile  systems.  Expertise  in  the  latest  scientific  devel- 
opments is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 

Aerophysus  laboratory;  Launch  and  reentry  aerodynamics,  heat  trans- 
fer. reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high -power  gas  lasers. 

Chemistry  and  Physics  Laboratory:  Atmospheric  reactions  and  atmos- 
pheric optics,  chemical  reactions  in  polluted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo- 
sensitive materials  and  sensors,  high  precision  laser  ranging,  and  the  appli- 
cation of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronic  s Research  Laboratory:  Electromagnetic  theory,  devices,  and 
propagation  phenomena.  Including  plasma  electromagnetics;  quantum  electronics, 
lasers,  and  electro-optics;  communication  sciences,  applied  electronics,  semi- 
conducting. superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging;  atmospheric  pollution;  millimeter  wave  and  far-infrared  technology. 

Materials  Sciences  Laboratory:  Development  of  new  materials;  metal 
matrix  composites  and  new  lormi  of  carbon;  test  and  exaluation  of  graphite 
and  ceramics  in  reentry,  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  environment,  application  of  fracture  mechanics  to  stress  cor- 
rosion and  fatigue-induced  fractures  in  structural  metals. 

S£ace  Sciences  Laboratory:  Atmospheric  and  ionospheric  physics,  radia- 
tion from  the  atmosphere,  density  and  composition  of  the  atmosphere,  aurorae 
and  airglow;  magnetospheric  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  the  magnetosphere;  solar  physics,  studies  of  solar  magnetic 
fields;  space  astronomy,  x-ray  astronomy;  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphere;  the  effects  of  optical,  electromagnetic,  and  particulate  radia- 
tions in  space  on  space  systems. 
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